Four-parameter families of complex Hadamard matrices of order six 
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In this paper we provide a general method to construct four-parameter families of complex 
Hadamard matrices of order six. Our approach is to write a 6-dimensional matrix as composed 
of four blocks, each one in the form of a circulant matrix. Use of Sylvester method of inverse 
orthogonal matrices leads to two constraints on the six parameters, and their fulfilment generates 
complex Hadamard matrices when the free parameters take values on the unit circle. We hope that 
the problem of mutually unbiased bases could be solved for six-level systems by using our results. 
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INTRODUCTION 



Construction of complex 6x6 Hadamard matrices has generated a considerable interest in the last years, this case 
being quite accessible, leading to a few numerical matrices and to several matrices depending on few arbitrary phases, 
see Refs. 0; 03 j 0~ 15 1 , [lj]]- 18 1, [H|. It was conjectured in paper Q that the above results are particular 



cases of a more general four-parameter family yet to be discovered. The number four is generated by the so called 
defect of a unitary matrix, see [2fj| |. 

It seems that until now there is only one attempt to construct such a family, see paper [l8j , where his author starts 
. with a dephased 6-dimensional matrix composed of four blocks and imposes some orthogonality conditions to get a 
O-f Hadamard matrix. However no explicit form for such a matrix is provided. 

The orthogonality concept is essential for Hadamard matrices and in the following we will use the method of inverse 
orthogonal matrices introduced by Sylvester in paper [l9j in the most general form. In this paper we will employ only 
a particular class of inverse orthogonal matrices, O — (oy ), namely those matrices whose inverse is given by 

o- 1 = (i/ 0lJ y = (i/ 0ji ) (i) 



' ^ \ where t means transpose, and their entries =^ o^- € C satisfy the relation 

OCT 1 = nJ n (2) 

When the Ojj entries take values on the unit circle, O -1 coincides with the Hermitean conjugate O* of O, and in this 
. ■ case relation ((2|) is the definition of complex Hadamard matrices. Accordingly the problem is to construct inverse 
C \ orthogonal matrices that will provide the corresponding complex Hadamard matrices. 

Our starting point is to write an arbitrary matrix as a four block matrix of the simplest form, as 



M 



A B 
C D 



(3) 



similar to that from paper [l9(, where A and B are circulant matrices of arbitrary dimension n. As concerns the C 
and D matrices they are given by C = B^ 1 and D = — A^ 1 , according to relation ([1}. The real problem is to find 
those constraints on A and B entries whose fulfilment transforms M into a complex Hadamard matrix, i.e. the new 
M matrix has unimodular entries and satisfies the relation 

MM' 1 = 2nl 2n (4) 

An important problem for complex Hadamard matrices is that of their equivalence. Nowadays one makes use of two 
different equivalence methods: the standard method and the unitary equivalence. The first one is currently applyed 
to real Hadamard matrices and has its origin in Sylvester paper who introduced the so called standard form for real 
matrices, whose entries of the first row and column are equal to 1. It was also extended to complex Hadamards. 

Its nowadays form says that two Hadamard matrices Hi and H 2 are equivalent, written as Hi = H 2 , if there exist 
two diagonal unitary matrices D\ and D 2 , and permutation matrices P\ and P 2 , such that 

H x = D X P X H 2 P 2 D 2 (5) 
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However the complex Hadamard matrices naturally belong to the class of normal matrices. A matrix N is normal 
if it commutes with its adjoint N* , i.e. it satisfies the relation N N* — N* N. For this class of operators the unitary 
equivalence takes a simple form and says that every normal matrix is similar to a diagonal matrix D, which means 
that there exists a unitary matrix U, such that 



N = UDU* 



(6) 



see [H p. 357, or [TTjj . 

Two important classes of normal operators are the unitary and self-adjoint matrices, such that for both these classes 
two matrices M\ and Mi are unitary equivalent iff they have the same s pec trum, or equivalently, the characteristic 
polynomials are the same up to a multiplicative constant factor, see [IF 



det(x I n — Mi/^/n) = det(x I n 



s pec trum, 
. HH, i.e. 



(7) 



where det is the determinant of the corresponding matrix. 

For unitary matrices the entries of the diagonal matrix D are unimodular, and for self-adjoint they are real numbers. 

To understand better the difference between the two equivalence methods we give a simple example of a selfadjoint 
Hadamard matrix whose entries from the first row and column are equal to 1. 
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The unitary equivalence implies the knowledge of its spectrum, which in this case is given by 

Sp(d 6 )= [-1 3 ,! 3 ] 



(9) 



where power means eigenvalue multiplicity. The above result shows that d$ is a complex Hadamard matrix, and in 
the same time a selfadoint matrix, its spectrum being unimodular and real. 
If we make use of the standard equivalence ([3]) the following matrix 



dei = 

is equivalent to ([5]), but its spectrum is given by 

Sp(d 6 i) = 
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(10) 



-l 2 l 2 



- V2 i + y/2 



V3 



V3 



(11) 



Thus when one makes use of the usual equivalence, see relation ([5]), the classical quantum mechanics will be into a 
big danger, because the dgi eigenvalues are unimodular, but not real. 

In the first case d§ = <ig, and in the second case dgi ^ ^61- Thus the change of rows and/or colums between 
themselves modifies the matrix symmetry, and accordingly its spectrum. 

The paper is organized as follows. In Sec. 2 we treat the simplest case, namely that when A and B matrices 
are 2-dimensional circulant matrices to see how the method works. In this simple case we show that the Sylvester 
orthogonality is equivalent to one condition on the four arbitrary parameters entering A and B matrices which leads 
to four non dephased and inequivalent complex Hadamard matrices. Similar results for n = 6 are given in Sec. 3. 
In Sec. 4 we show how the results obtained in sections 2 and 3 lead to higher order Hadamard matrices. The paper 
ends by a few conclusions. 
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2. FOUR-DIMENSIONAL HADAMARD MATRICES 



The current opinion about four-dimensional Hadamard matrices is that up to equivalence there is only one matrix, 
namely that found by Hadamard Q , which is written in a dephased form. As we said this approach has its source in 
Sylvester paper who introduced the so called standard form for real ±1 matrices. This equivalence was essentially 
extended to all complex Hadamard matrices. But if we look at Hadamard matrices as unitary matrices, then their 
equivalence is given by their spectra, according to papers (l6| . or [TTj . and the problem changes significantly. 

If one takes M of the form 



Ma = 



a 
-b 

_L 

c 
_L 

d 



b C 

-d 
1 



a 



the relation @ leads to the following orthogonality condition 

(be + ad)(ac — bd) 







(12) 



(13) 



and if it is satisfied, one gets an orthogonal matrix. As it is easily seen the above equation has eight solutions, e.g. 
a = b = — — , c = — , etc, but the unitary equivalence principle restricts them to four nonequivalcnt ones. 
If we start with a = — the orthogonal matrix has the form 



Hn(b, c, d) — 



b 

bd 

c 

d 



bd 

_ 1 
b 



bd 



(14) 



When the parameters b, c, d € T 3 , the three dimensional torus, the matrix (fT4|) is complex Hadamard and its spectrum 
is given by the solutions of the equation 



eci(6, c,d;x) = x 4 



(b 2 d 2 -c 2 ) 3 (c 2 +d 2 )(c 2 +6 4 d 2 ) 



9c 2 d 2 



(b 2 d 2 



bed 



4b 2 c 2 d 2 



bed 



1 = 



(15) 



which depends on three independent phases. The above equation can be easily solved. However since (|15[) is a 
reciprocal equation by the substitution 1/x — x = y it is reduced to a second order equation. This property is almost 
generic for any complex Hadamard matrix obtained by using relation such that its spectrum can be explicitly 
computed for dimensions n = 4, 6, 8. However there are few cases for dimensions n > 6 when only the numeric 
spectrum can be obtained, and they will be discussed in an other paper. 

Indeed the above substitution transforms (|15l) into the following simpler equation 



, u J x 2 (b 2 d 2 -e 2 ) (c 2 + d 2 )(e 2 +b 4 d 2 ) 
eci(6, c, d; x) = y 2 - — — -y + , n = 



bed 



Ab 2 c 2 d 2 



and its solutions are 



Via 



(b 2 d 2 -c 2 )± y/-(l + b 2 ) 2 C 2 d 2 



2bcd 

which provide the spectrum of equation (11 5[) 

b 2 d 2 - c 2 - y/-(l + b 2 ) 2 c 2 d 2 ± ^jl6b 2 c 2 d 2 + (c 2 -b 2 d 2 + ^-(l + tfflNpy 



Xl.2 



X 3 .i 



Abed 



b 2 d 2 - c 2 - J- 



v /-(l + 6 2 ) 2 c 2 d 2 ± J 16 b 2 c 2 d 2 + (c 2 - b 2 d 2 - y/-(l + 6 2 ) 2 c 2 "5y 



Abed 



(16) 



(17) 



(18) 
(19) 
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For example the parameters b = 1, c = i, d = —i provide the spectrum 



Sp(ff 4 (l, *,-*)) 



Zl,2 



2^3,4 



± V3 



where in the following i = y/—l. 

Equation (JPojl for 6 = ±1, c= e T , (i = eT transforms into 



2x 4 + 2iv / 2a; 3 - 3a; 2 - 2iV2x + 2 = 



and its spectrum is given by 

Sp(H 4 (l,e^,e^)) = 
The standard form of H 4 matrix is 



(2- y/2)i± J2(5 + 2y/2) -(2 + y/2)i± J2(5 - 2y/2) 



H 4a = 



1111 

1 -q q -1 
1-1-1 1 
1 q -q -1 



where q — dr/c and its spectrum is given by 



Sp(H 4a (q)) = 

with g an arbitrary phase. From the above formula one hnds 



1 + g + \/j - 14g + g 2 " 1 + g - - 14g + . 

' ' 4 ' 4 



Sp(fl4,(l)) = 



Sp(H 4a (i)) 



-M,- 



l + iV3 -1 + «V3 



1,1, 1 + *_ L ( 1 -0>/7 



5p(^4a(-l)) = [-1 2 ,! 2 ] 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



The spectrum of the original Hadamard matrix found in paper Q is obtained from formula (|24l) by the change q — > —q. 
If one interchange the third and the fourth rows of (|23|) between themselves the spectrum is given by [—1, 1, 1, — q]. 
This example shows that by bringing a unitary matrix to its dephased form, or by permuting rows and columns 
between themselves, the obtained matrices are in general not unitary equivalent with the matrix one started with. 

The parameter choices, c— ^,d— ^,0 — — and d = — — , lead to only one matrix, H 4 2(a, 6), which depends 
on two parameters, and its spectral equation is 



n 2 - 1 I- 2 
ec2(a, 0; x) = x x 



{a 2 + b 2 ){l + a 2 b 2 ) - 8a 2 b 2 



a 2 - 1 



-x + 1 =0 



a Aa 2 b 2 a 

The matrix obtained with the condition b = 2£ H 4 ^(a, c, d), has the spectrum given by the equation 



ec3 (a, c, d; a;) = x 4 



1 3 (c 2 + d 2 )( a 4 c 2 + d 2 )-8a 2 c 2 rf 



4a 2 c 2 <i 2 



2 - 2 ^ 2 2 a 2 -l 

2 1 x + 1 = 



■a; 



The choice a = — ^ leads to H 44 (a, c, d) matrix with the spectral equation 

ec 4 (b, c, <i; ar) — x i + 
The choice & ; 



6 2 c 2 -d 2 , (c 2 + d 2 )(6 4 c 2 + d 2 )- 86 2 c 2 d 2 2 6 2 c 2 -d 2 

1 X 1 ; ; X + 1=0 



-X 



bed 4b 2 c 2 d 2 
— provides the matrix H 4 $(a, c, d) whose spectral equation is 



bed 



(26) 



(27) 



(28) 



5 



t-. J , 4 {a 2 -I) 3 (c 2 + d 2 )(c 2 + a 4 d 2 )-86 2 c 2 d 2 2 (a 2 - 1) 1 n 

ec 5 (b,c,d;x) = x 4 - -x 3 + ^ A , 9 „ ' x 2 + ± '-x + 1 = 29) 

a 4a^d^ a 

equation which is similar to equation (|27p . 

A more careful analysis shows that the spectra of equations (|2"o| and (|2"T)) coincide on the surface d 2 = a 2 b 2 c 2 , i.e. 
the matrix (a, b) is a particular case of the matrix H^{a, c, d). 



SIX-DIMENSIONAL HADAMARD MATRICES 



In order to see a few subtleties of case n = 6 we start with the Bjorck and Froberg 6-dimensional matrix, (BJ , written 
in the form provided by Haagerup, Q 



BF = 
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(30) 



The orthogonality relation (j4)) leads to one constraint 



whose solutions are 



2d A - 2d + 1 = 



_ 1 ± iy/2 3 1 / 4 - V3 1 ± 3 1 / 4 + V3 

"1.2 — n 1 "3,4 — 7: 



(31) 



(32) 



The first two solutions are complex and unimodular and lead to one complex Hadamard matrix. The d3.4 solutions 
being real numbers different from ±1 lead to two different orthogonal matrices. The spectral equation is the same for 
both d\ 2, and its form is 



c 6 -V6x 5 +3x 4 -2V2x 3 + 3x 2 -V6x + 1 =0 



(33) 



and the eigenvalues are given by 



Sp(BF) 



l±i VQ-V2 + 23^ 1 - 3 V&-V2 + 2Si / 1 - ^2 V3 - 3 

VT i ±Z V 2^ J± i l + 2 



(34) 



The entries of the two Sylvester orthogonal matrices are not unimodular so they do not generate Hadamard matrices. 
The dephased form of the matrix is 



bf = 
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1 
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1 
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-d 2 
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1 

d' 2 


1 

d 2 
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-d 3 
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-d 2 


-d 


-1 



(35) 



and if we substitute in it, e.g. the di— root, its spectrum is Sp(bf) = [— l 3 , l 3 ] . Thus the matrices (|30[) and (|35[) arc 
not unitary equivalent. 
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The 6-dimensional matrix is of the form ([3]) whith A and B blocks 3-dimensional circulant matrices, and explicitly 
it is written as 



a 


b 
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d 




f 1 
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e 
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a - 



(36) 



The orthogonality relation (Q| leads to two constraints 



abcd 2 e + a 2 bdef + b 2 cdef + ac 2 def + abce 2 f + abcdf 2 =0 
abcde 2 +abcd 2 f + ab 2 def + a 2 cdef + bc 2 def + abcef 2 = 



(37) 
(38) 



and when they are satisfied the matrix (|36[) gets orthogonal. Each equation is quadratic in a, 6, c, <i, e, / parameters 
such that one could solve one of them with respect to a chosen parameter. For example we can use equation (|37[) to 
solve the constraint on / parameter, that provides two solutions for it 



f± 



-e(a 2 bd — b 2 cd — ac 2 d — abce) ± yj —Aa 2 b 2 c 2 d 3 e + e 2 (a 2 bd + b 2 cd + ac 2 e + abce) 2 

2abcd 



(39) 



where the subscript of f± is related to the ± sign which appears in front of the square root in relation (I39|) . When 
both solutions f± are introduced in the equation (|38|) one gets only one condition 

- a b c d 3 - a b 2 d 2 e - a 2 c d 2 e - b c 2 d 2 e + a 2 b d e 2 + b 2 c d e 2 + a c 2 d e 2 + a b c e 3 = (40) 

If one solves equation (j38]) with respect to / parameter, and if these solutions are introduced in relation (1371) one 
gets again the relation (|4T7]). 

This equation has a quadratic dependence in a, &andc parameters, and a cubic dependence in dande parameters. 
Thus we can solve it with respect to any of the first three parameters by getting six solutions for a±, b±, and c±, and 
six solutions, di,2,3 and ei.2,3, for cubic parameters. For example a± have the form 



a± 



(cd + be)(bd 2 -ce 2 )± y/{cd + be) 2 {bd 2 - ce 2 ) 2 ~ iabcd 2 e 2 (cd - be) 



2de(be — cd) 

If we solve equation (|37[) with respect to parameter a one finds 



a± 



-c(bd 2 e + cdef + be 2 f + bdf 2 ) ± v /-4 6 3 cd 2 e 2 f 2 + c 2 (bd 2 e + cde f + be 2 f + b d f 2 



When a + is introduced into equation 



2abcd 

one gets the following condition 



- be 2 d 2 e + b 2 cde 2 + b 2 cd 2 f + b 3 de f - c 3 de f - be 2 e 2 f - be 2 d f 2 + b 2 ce f 2 = 
and if we solve it for / parameter, its solutions are 



f± = 



(cd + be)(bd 2 - ce 2 ) ± y/ {cd + be) 2 (bd 2 - ce 2 ) 2 - Aabcd 2 e 2 (cd - be) 

2de(be - cd) 



(41) 



(42) 



(43) 



(44) 



and two solutions for d± and e±, and three solutions for ^1,2,3 and ci^.3- 

The a± solutions from the first case, (|4H) . coincide with the f± solutions from the second case, (|4^|) . whilst the 
solutions for f± from the first case, (13^1) . are different from the solutions for a±, (|4"2"]) . Thus the above coincidence 
suggests that there are only four independent matrices instead of eight for each pair (a, /), (d, e), etc. To prove it one 
has to find the spectral equation for each case, which is not a simple problem. 

In contradistinction with the precedent four dimensional case the obtained matrices are large having a complicated 
dependence on a± and f±, as it is seen from relations (|3"9"|) . (|4"Tj) . (|4"2")l . (|4"4")l . such that they cannot be written on an 
usual paper sheet. 
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To have an idea about their spectra we make use of numerical values for the remaining parameters. For example if 
one chooses for free parameters the following values, b = 1, c = i, d = e~ , e = — 1, one get four numerical Hadamard 
matrices, and their spectra are given by the equations 



ec(f+,a+,x) = x 6 +iV3x 5 - 3x 4 - ~iy86 + ^Lll x 3 + 3x 2 + iV3x - 1 = (45) 



86 + 


32\/2 
3 


86- 


32%/2 
3 


86 + 


32a/2 
3 



ec(f+,a-,x) = x 6 -iV3x 5 -Sx 4 + ^i\l86-^—x 3 + 3x 2 -iV3x-l = (46) 



(47) 



ec(/_,a_,x) = x 6 - iV3x 5 - 3 x 4 + ~i<t /86 - ^l1 x 3 + 3 x 2 - iVix - 1 = (48) 

O y O 

where f± and a± are provided by relations and respectively (|4ip. By the substitution y = i — x the above 

equations transform into reduced equations 



^ 9 -18+ J2 129+16V2 ^ 9 -18 + 2J2 129-16^2 

y 3 -iV3y 2 + i V . = Q, y 3 + iV3y 2 - t i = =0 (49) 

V3V3 3V3 

If one makes use of the second set (o±, f±) one gets 



ec(a+,f + ,x) = x 6 + iV3x 5 - 3x 4 - -\l 86 + x 3 + 3x 2 + iV3x - 1 = (50) 



ec(o_,/+,x) = a; 6 - i%/3x 5 -3x 4 + ^1/86- ^^x 3 + 3x 2 -iv^x- 1 = (52) 



86 + 


32^ 
3 


86- 


32V^ 


3 


86- 


32^2 


3 



ec(o+,/_,x) = x 6 -i%/3x 5 -3x 4 + ^i/86 ^— x 3 + 3x 2 - iV3x - 1 = (51) 



ec(a-,f-,x) = x 6 + iV3x 5 - 3x 4 - ^/86+ 5^x 3 + 3x 2 + iV3x - 1 = (53) 



and the reduced equations are the same as in the preceding case. 

The above numerical examples could suggest that the pairs (/+,a±) and (a_,/±) lead to different spectra. Taking 
into account that the relations (I57|) and (|38p can be solved for all pairs of parameters it follows that, in general, it 
could be 4 x Cg = 4 x 15 = 60 inequivalent matrices. 

Supplementary solutions can be obtained from equations similar to (|40p by solving them with respect to cubic 
parameters; in the case of equation (|40|) . the d or e parameters. In all this cases the solutions are more complicated, 
and the corresponding matrices also. For example if we solve equation (140)) with respect to e parameter one expects to 
find three solutions. However in some numerical cases the spectral equation does not depend on f±, or 61,2,3 solutions 
of the cubic equation. For example the parameters, a = 1, b = e 2 ^, c = e^ 1 , d = 1, lead to the spectral equation 
whose spectrum is simple 



:(/±,ei, 2 , 3 ) = (x 2 ~ 1)(^ 4 + 1) = 0, Sp(/ ±) ei, 2 , 3 ) = 



1 + i 1-i 



(54) 



If we solve equation (1401) with respect to d parameter on finds for a = 1, b = e 2 " , c = e 4 3* , e = 1 the following 
spectrum 

W±,di,a,3) = [-1 3 , I 3 ] (55) 



which coincides with the spectrum of the matrix (|35[) . Thus the dephased form of Bjorck and Froberg matrix is 
unitary equivalent to a particular case of the matrix (|36[) when a second parameter is obtained from a cubic equation. 
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Thus for each parameter, a, b, c, d, e, /, obtained from the solutions of cubic equations, one gets at least a single 
Hadamard matrix, such that the number of all nonequivalent matrices could be at least 60 + 6 = 66. 

The above matrices are a rich source for three-parameter Hadamard matrices. The simplest case will be to equal 
two parameters; however the obtained matrix is also large. Many simpler matrices which depend on three parameters 
do exist and in the following we give a few examples. 

Looking at relation (|4Tj) one easily sees that if the following equation is satisfied 



(be + cd)(bd 2 - ce 2 ) = 



(56) 



we can obtain three-parameter Hadamard matrices as follows. With the solution b = —cd/e from the first factor (J5j 
one gets the Hadamard matrix 



£>6i(c,d,e) 



\/c 4 d 5 < 
c d 2 e 



c d 



Vc 4 d 5 e 
cd 2 e 



Vc 4 d 5 e 
cd' 2 e 



c V c 4 d 5 e 

e 



cy_ c 4 d 5 e 



cV c 4 d 5 e 



cVc 4 d 5 < 



Vc r dF< 

e 

c d 
1 



cV c 4 d 5 e 

_ 1 

c 

cd 2 e 



cVc 4 d 5 < 

e 
d 

C d 
1 



Vc 4 d s ; 

e 

c d 



Vc 4 d 5 e 



which comes from the pair (/ + , a + ) in this order. Its reduced polynomial has the form 

^/2{c 2 d-e)V^d^ 2 (c 4 d 2 + e 2 )y 



Pol(f+, a+)=y - 



c 2 de 



2 d(d 3 - e 3 )(c 2 d + e) 3 + 4(c 6 d 3 - e 3 )vW^ 



6 c 5 d 4 e 3 



The next matrix comes from the pair (/+, a_) 



D 62 (c, d, e) 



\/e 4 d 5 e 
cd 2 e 

C 

c d 



Vc 4 d 5 e 
cd 2 e 



Vc 4 ri 5 6 
c d 2 e 



cv c 4 d 5 e 



cV c 4 d 5 e 



cVc 4 d 5 e 



cV c 4 d 5 e 



Vc 4 d 5 e 
e 

c d 
1 



V- 



cVc 4 d 5 « 

_ 1 

c 

c d 2 e 



Vc 4 d 5 . 
e 

e d 



cV c 4 d 5 e 

e 
d 

e 

c d 

_ 1 

c 

erf 2 e 
\/c 4 d 5 e 



and its reduced polynomial is 



(57) 



(58) 



(59) 



pol(f+, o_) = y 3 + 



^(^d - ejVc 4 ^ 2 (c 4 d 2 +e 2 )y 



c 3 d 3 e 



c 2 de 



c 2 d(d 3 - e 3 )(c 2 d + e) 3 - 4(c 6 d 3 - e 3 )yWt 
6 V6"c 5 d 4 e 3 



(60) 



From the reduced equations (1551) and one sees that the above two matrices are not unitary equivalent. 

There are two more matrices generated by the pairs (a+, /+) and (a + , /_), and for them we write only their reduced 
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polynomials 



pol(a+, /+) = y 3 + 



/3/2(c 2 d-e)V^dP^~- 



c 6 d-> 



(c 4 d 2 



c 6 d 6 y c 2 de 

- 5 d 4 (d 3 - e 3 )(c 2 d + e) 3 - 4(c 6 d 3 - e 3 )^/ 3 ^/^ 
6V6c 8 d 7 e 3 



(61) 



pol(a+, f- 



^J2{c 2 d - e)^¥^f 2 (c 4 d 2 + p2) 



3 rF> 



c 2 de 



~ 5 d 4 {d 3 - e 3 )(c 2 d + e) 3 + 4(c 6 d 4 - e 3 )V^¥^-^f 
6^^ d 7 e 3 



(62) 



The reduced four polynomials, (|58[) and (|60p - (|62[) . are a strong argument that the four-parameter matrices have 
only four independent solutions for each pair a± and f±, and those similar to. The above three-parameter matrices 
obtained from four-parameter familes show that the number of these matrices is at least 60. 

From the second factor of relation (f56"|) one get similar results. For example the solutions b — ce 2 /d 2 and c = bd 2 /e 2 
lead also to two matrices which are a litlle more complicated than the D 61 and D 62 matrices. The other solutions 
d = ±y/ ce 2 Jb and e = ±\/b d 2 / c could lead in principle to four matrices. Their spectral equations are also more 
complicated. 

All these matrices are different from the other three-parameter matrices found in literature, e.g., see paper [liij . 
They cannot be "simplyfied" to simpler forms; in some cases the analytic continuation can be done by "hand" , but 
the spectra of the new matrices are in general different from the spectra of matrices one started with. From the above 
matrices one get many matrices depending on two parameters by making equal two of them. 

As concerns the number of three-parameter matrices the second factor of relation (|56p generates also four different 
matrices similar to matrices ([5T)) - (f6"2")) such that the total number of such matrices is at least 8 x 15 = 120. 

Let see what happens when the matrix Mq is brought to its standard form. One gets 



Mb, = 



1 


1 


l 


l 


1 


2 

a 
be 


a b 

c 2 


El 

cd 


atl 

ce 


ae 
cf 


ac 
b 2 


a 2 
be 


ac 
bd 


af 
be 


ad 
bf 


ad 
bf 


ad 

ce 


-1 


ad 

ce 


ad 
bf 


ac 
bd 


ac 
cf 


ae 
bd 


-i 


ae 
cf 


a£ 

be 


nl 

cd 


_«£ 

cd 


af 
be 


-1 



(63) 



The constraints satisfied by Mq s matrix are the same as for the matrix Mq, i.e. the equations (pTTf and (|38l) . and 
after their fulfilment its spectral equations have to be in general different from those generated by the Mq matrix. 

As an example we generated matrices by using the pairs f± and a± into the matrix Dq s , and we made two choices, 
b = —cd/e and b = ce 2 /d 2 , to obtain matrices similar to matrices whose reduced polynomials are (|58|) and ([6"0l) - ()62[) . 
The unexpected result was that there is only one polynomial, and one eigenvalue set 



ec(M, 



(is) 



Vr 5 



1=0, Sp(M es 



-1,1, 



1 + 



(64) 



which does not depend on f± and a±, showing that the dephased and its undephased forms are not unitary equivalent. 

An explanation could be the following. The coefficient of the x 3 power in the spectral polynomial of the matrix 
which depends on the four free parameters, 6, c, d, e, has as factor the relation (|56[) multliplied by a huge expression. 
This explains the absence of the x 3 power in equation (|64[) . but does no explain the other constant terms. 
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4. HIGHER DIMENSIONAL HADAMARD MATRICES 



For n = 8 the analog of matrix (|36|) is the following 



M 8 = 



a 


6 


C 


d 


e 


/ 


.9 




d 


a 


6 


c 


ft 


e 


./' 


g 


c 


d 


a 


b 


9 


/i 


e 


f 


b 


c 


d 


a 


f 


.9 


h 


e 


i 


i 


i 


i 


l 


l 


i 


i 


e 




9 


/ 


a 


d 


c 


b 


1 


_L 


1 


_L 


1 




1 


_ 1 


/ 


e 


h 


.<; 


6 


a 


d 


c 


1 


1 


1 


l 


1 


1 


1 


1 


9 


/ 


e 


h 


c 


6 


a 


d 


1 


1 


1 


1 


1 


1 


1 


_ J_ 


ft 


9 


f 


e 


d 


c 


6 


a _ 



and relation (j4]) leads to three constraints 

abcde 2 fg + a 2 bcefgh + b 2 cdefgh + ac 2 defgh + abd 2 efgh + abcdf 2 gh + abcdeg 2 h + abcdefh 2 = (66) 
abode f 2 g + abode 2 fh + ab 2 cefgh + a 2 bdefgh + bc 2 defgh + acd 2 efgh + abcdfg 2 h + abcdegh 2 — (67) 
abcdefg 2 + abcdef 2 h + abcde 2 gh + abc 2 efgh + ab 2 defgh + a 2 cdefgh + bcd 2 efgh + abcdfgh 2 — (68) 

Solving equation ([6"6"]) with respect to /i parameter one get two solutions. When these solutions are introduced in 
equation (|68p one gets one equation which does not depend on the square root, similar to the case n = 6 

(a 6ce/ + 6 cc?e/ + ac def + abd ef + abcdf + abcdeg)g — 

(abcdf 2 + abcdeg + abc 2 fg + ab 2 dfg + a 2 cdfg + bed 2 fg)e 2 (69) 

If the same solution is introduced into equation (|6T[) one gets an equation that depends on the square root entering 
the h parameter. By multiplying it with its conjugate one gets a huge polynomial, P, that depends on powers, 
a 4 , 6 4 , ... , g 4 , in the corresponding monomials. If one solution of the equation (|69| is introduced into the above P 
polynomial one finds an equation that depends on square root and the third power of the square root expressions. 
Even if such an equation could be transformed into a polynomial, each power of the remaining parameters will be 
greater than eight, such that one cannot solve analytically such an equation. 

On could say that the form of (|65p matrix is not appropiate for getting Hadamard matrices; we tried a few other 
forms but the obtained constraints are at least four, such that in the best case the resulting matrix will depend on 
the same number of parameters as matrices obtained from ([36]). 

Thus in this case we have to make use of the alternative formula, see equation (8), of the form 



where A and B are Hadamard matrices, and D is a n-dimensional diagonal matrix containing phases. 
By using the results from section 2 we can take two 4-dimensional matrices 





' bd 

c 


b 


C 


d~ 




f 


9 


h 


fh - 

9 




-b 


bd 


-d 


c 




-g 


f 




h 


A 4 = 




c 








9 






i 


1 


c 


i 


i 


2_ 


1 


1 




c 


d 


bd 


b 




h 


fh 


/ 


a 




1 


1 


1 


c 






1 


_1 


l 




. d 


c 


b 


bd . 




L fh 


h 


a 


f J 



(71) 



and construct with them an eight-dimensional matrix by using relation (|70|). It has the form 
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D 8a (b, c,d, f,g,h) 



~ bd 

c 


b 


C 


d 


f 


9 


h 


fh - 

3 


-b 


bd 

c 


-d 


c 


-.9 


f 


a 


h 


i 

c 


1 
d 


c 

bd 


i 

b 


l 

h 


g 

fh 


i 

/ 


_L 

9 


1 

d 


1 

c 


1 

h 


c 

bd 


9 

fh 


l 

h 


_ i 

9 


1 
/ 


bd 

c 


b 


C 


d 


-f 


~9 


-h 


_£h 

9 


-b 


bd 

c 


-d 


c 


9 


-f 


III 

9 


-h 


1 

c 


1 
d 


c 

bd 


i 

b 


l 

h 


y 

fh 


1 

f 


_i 

9 


1 

. d 


1 

c 


1 

b 


c 

bd 


s 

fh 


1 

h 


1 

a 


1 

/ J 



(72) 



The above matrix depends on 6 free parameters and if the i?-matrix is multiplied by D the new matrix depends on 
10 parameters. 

The spectral equation of the matrix (|72[1 can be written but it is too long, so we chose some particular values for 
all the six parameters, 6=1, c = i ,d = —i, f = g = h = — 1, which generate the matrix 



£> 8 i(l,i,-i,e",e 



The reduced equation has the form 



-1) 



y -iy + 



-1 


1 


i —i 


-1 


-1 


i i 


-1 


—i 


1 1 


i 


—i 


-1 1 


-1 


1 


i —i 


-1 


-1 


i i 


—i 


—i 


1 1 


i 


—i 


-1 1 


V2 


.6 

— ; 


+ V2 



l-H 


l-i 


V2 


V2 


l-i 


l+i 


V2 


V2 


-1 


—i 


i 


-1 


l+i 


l-i 


V2 


V2 


l-i 


l+i 


V2 


V2 


1 


i 


—i 


1 



1 + V2 



-1 


—i 


i 


-1 


l-i 


l+i 


V2 


V2 


l+i 


l-i 


V2 


V2 


1 


i 


—i 


1 


l-i 


l+i 


V2 


V2 


l+i 


l-i 


V2 


V2 J 



(73) 



= 



(74) 



whose solutions are 



Vi 



iy/2, 



>J2 



V2 



.{V2 + VTo) 



?74 



.(V2- VTo) 



(75) 



and the spectrum of matrix (1731) is given by 



Sp(D S i{l,i,-i,e" ,e " ,-1) 



i± 1 

^1,2 = ^3,4 = 



(2 + a/2) ± Vl0-4\/2 



V2 



±2^13 - -y/5 - «\/2(l + V5) i(\/T0- ^±2^/13 + V5 
X 5,6 = s , ^7,8 = 3 



(76) 



One problem is to find how many 8-dimensional matrices, similar to (1721 can be found. In section 2 we have shown 
that the relation (| 13[) has eight solutions. Thus all these matrices can be used as matrix A, and by changing the 
parameters as in the above matrix B one arrives at 64 matrices similar to (|72[). 

The relation ([70]) can be used again with the matrix (1721) to obtain matrices of dimensions i?2™8 for n = 1, 2, . . . ; 
for example the matrix Hie will depend on 10+10+8=28 free parameters, etc. 

The same thing can be done with the 6-dimensional Hadamard matrices. All the 4-parameter six-dimensional 
matrices can be used to obtain 12-dimensional Hadamard matrices which depend on fourteen parameters. These 
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matrices are also large and cannot be written down, but numerical matrices can be. The spectral equation will be 
12-dimensional and its reduced equation is 6-dimensional, so one cannot find for it analytic solutions. A numerical 
example is the following 



A. 



V2 

i 
1 

i 

-1 
l+i 

V2 



1 


i 


—i 


-1 


1-i 
s/2 


1 


l-i 
V2 


—i 


i 


l-i 
V2 


-1 


l-i 


l+i 
V2 


-1 


l+i 
V2 


i 


i 


l+i 


-1 


l+i 
V2 


-1 


i 


i 


-1 



s/2 
-1 

— % 

-1 



l+i 
V2 



Bq = 



2 1 / 4 

l+i 
V2 

yr+i 

2 l/4 

i+i 
-l 



i+i 

1 

t-y/l+i 

w 
-1 

yr+i 

21/4 

i+i 

72 



2 1 / 4 

l+i 

V2 

1 

l+i 

V2 

-1 

VT+i 

W 1 



VT+i 

-1 

l-i 
V2 

-1 

l-i 

3/2 



(1 + 



%/2 

VT+i 
-1 

(l+i) 3/2 
2 3 / 4 

-1 

l-i 

71 



-1 

l-i 

2 l/4 



l-i 
V2 

(l+i) 3/2 
2 3 / 4 

-1 



(77) 



The spectrum of the corresponding numerical D\i matrix generated by Aq and Bq can be found numerically by 
using, e.g. Mathematica, and in this numerical case the spectrum is simple. As in the preceding n = 8 case, for 
ii > 12 one can use relation (|70p to obtain Hadamard matrices with dimensions H<x n ii for any integer n > 1. 



5. CONCLUSION 



Starting with a carefully chosen form of the 6-dimensional matrix, (|36p . we got many 6-dimensional Hadamard 
matrices that depend on four parameters. The peculiarity of all of them is that they are long matrices which can be 
stored and viewed only in an electronic format. In fact after the substitution of two chosen parameters, e.g. b and 
/, the complex Hadamard matrix Mq — Mq(cl, c, d, e) depends on four independent parameters. However many of its 
entries become so long, that from a typographical point of view Mq can hardly be presented as a matrix, but it can 
be used in an algebraic computer program. 

The number of four- and three-parameter matrices being quite big there is a hope that the unsolved problem of 
mutually unbiased bases could be solved for six-level systems. This problem received a special attention of many 
groups, but despite substantial efforts until now no relevant breakthrough was provided, although papers appear 
continuosly, the last being [l(| and Q. 

Looking at the form of Mq matrix one sees that by our procedure one cannot obtain all the the 6-dimensional 
matrices that depend on four parameters, one of them is the Agaian matrix, see [l|, page 112. Thus other approaches 
are necessary to find all those matrices which include Agaian matrix. 
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